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1 

Qh A research program within the scope of theories on "Emergent Quantum Mechanics" is presented, 



which has gained some momentum in recent years. Via the modeling of a quantum system as a 
non-equilibrium steady-state maintained by a permanent throughput of energy from the zero-point 



rS^ vacuum, the quantum is considered as an emergent system. We implement a specific "bouncer- 

+L walker" model in the context of an assumed sub-quantum statistical physics, in analogy to the 

s 

g results of experiments by Couder's group on a classical wave-particle duality. We can thus give 

an explanation of various quantum mechanical features and results on the basis of a "21 st century 

classical physics", such as the appearance of Planck's constant, the Schrodinger equation, etc. An 

> 

Q\ essential result is given by the proof that averaged particle trajectories' behaviors correspond to a 

specific type of anomalous diffusion termed "ballistic" diffusion on a sub-quantum level. 
^- It is further demonstrated both analytically and with the aid of computer simulations that our 

o 

model provides explanations for various quantum effects such as double-slit or n-slit interference. 
We show the averaged trajectories emerging from our model to be identical to Bohmian trajectories, 



^ albeit without the need to invoke complex wave functions or any other quantum mechanical tool. 

H 

Finally, the model provides new insights into the origins of entanglement, and, in particular, into 

the phenomenon of a "systemic" nonlocality. 
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1. INTRODUCTION 

After many attempts to frame wave-particle duality on a single "basic" explanatory level, 
like on pure formalism (Heisenberg,. . . ), pure wave mechanics (deBroglie,. . . ), or pure par- 
ticle physics (Feynman,. . . ), for example, it is instructive to ask the following question: is 
perhaps the quantum a more complex dynamical phenomenon? One motivation for raising 
this question comes from the fact that we are presently witnessing a historical change with 
respect to wave-particle duality in double-slit interference. With regard to the latter, the 
standard attitude during practically the whole of the 20th century can be characterized by 
Richard Feynman's famous claim that it is "... a phenomenon which is impossible, abso- 
lutely impossible to explain in any classical way." In the present, however, we are facing the 
beautiful experiments by Yves Couder's group in Paris which do show wave-particle duality 
in a completely classical system, i.e., with the aid of bouncers/walkers on an oscillating fluid. 
(For a first impression, see the video with Yves Couder [J].) 

In fact, with our approach we have in a series of papers obtained from a "classical" 
dynamical model several essential elements of quantum theory [1-8]. They derive from the 
assumption that a "particle" of energy E = Hoo is actually an oscillator of angular frequency 
oj phase-locked with the zero-point oscillations of the surrounding environment, the latter 
of which containing both regular and fluctuating components and being constrained by the 
boundary conditions of the experimental setup via the buildup and maintenance of standing 
waves. The "particle" in this approach is an off-equilibrium steady-state maintained by the 
throughput of zero-point energy from its "vacuum" surroundings. In other words, in our 
model a quantum emerges from the synchronized dynamical coupling between a bouncer 
and its wave-like environment. This is in close analogy to the bouncing/walking droplets in 
the experiments of Couder's group [ >-12], which in many respects can serve as a classical 
prototype guiding our intuition. Our research program thus pertains to the scope of theories 
on "Emergent Quantum Mechanics" . (For the proceedings of a first international conference 
exclusively devoted to this topic, see Grossing (2012) [ ]. For their original models, see in 
particular the papers by Adler, Elze, Ord, Grossing etal, Cetto etal, Hiley, deGosson, 
Bacciagaluppi, Budiyono, Schuch, Faber, Hofer, 'tHooft, Khrennikov etal, Wetterich, and 
Burinskii.) We have recently applied our model to the case of interference at a double slit 
1 7], thereby obtaining the exact quantum mechanical probability distributions on a screen 



behind the double slit, the average particle trajectories (which because of the averaging are 
shown to be identical to the Bohmian ones), and the involved probability density currents. 

Moreover, already some decades ago, Aharonov etal. [14, 15] investigated the problem 
of double-slit interference "from a single particle perspective" . Their explanation is based 
on what they call dynamical nonlocality, and they therewith answer the question of how a 
particle going through one slit can "know" about the state of the other slit (i.e., being open or 
closed, for example). This type of "dynamical" nonlocality may lead to (causality preserving) 
changes in probability distributions and is thus distinguished from the kinematic nonlocality 
implicit in many quantum correlations, which, however, does not cause any changes in 
probability distributions. Chapter 3 of the present review extends the applicability also of 
our model to a such a "dynamical" scenario, which in our terminology rather relates to a 
"systemic" nonlocality, as shall be explicated below. 

This paper is organized as follows. In Chapter 2, some main results of our sub-quantum 
approach to quantum mechanics are presented. We begin by constructing from classical 
physics a model that explains the quantum mechanical dispersion of a Gaussian wave packet. 
As a consequence, we obtain the total velocity field for the current emerging from a Gaussian 
slit, which can easily be extended to a two-slit (and in fact to any n-slit) system, thereby 
also providing an explanation of the famous interference effects. With the thus obtained 
velocity field, one can by integration also obtain an action function that generalizes the or- 
dinary classical one to the case that is relevant for the reproduction of the quantum results. 
Effectively, this is achieved by the appearance of an additional kinetic energy term, which we 
associate with fluxes of "heat" in the vacuum, i.e., sub-quantum currents based on the exis- 
tence of the zero-point energy field. Finally, it is shown how the Schrodinger equation can 
be derived in a straightforward manner from the said generalized action function, or the cor- 
responding Lagrangian, respectively. In Chapter 3, then, we deal with the above-mentioned 
problem of interference from the single particle perspective, obtaining results which are in 
agreement with those presented by Aharonov etal. [14, 15], but differing in interpretation. 
Finally, in Chapter 4 computer simulations are presented to show some examples of how our 
sub-quantum approach reproduces results of ordinary quantum mechanics, but eventually 
may go beyond these. 



2. THE QUANTUM AS AN EMERGENT SYSTEM: A SUB-QUANTUM AP- 
PROACH TO QUANTUM MECHANICS 

We transfer an insight from the Couder experiments into our modeling of quantum sys- 
tems and assume that the waves are a space-filling phenomenon involving the whole ex- 
perimental setup. Thus, one can imagine a partial decoupling of the physics of waves and 
particles in that the latter still may be "guided" through said "landscape", but the for- 
mer may influence other regions of the "landscape" by providing specific phase information 
independently of the propagation of the particle. This is why a remote change in the ex- 
perimental setup, when mediated to the particle via de- and/or re-construction of standing 
waves, can amount to a nonlocal effect on a particle via the thus modified guiding landscape. 

In earlier work [■] we presented a model for the classical explanation of the quantum 
mechanical dispersion of a free Gaussian wave packet. In accordance with the classical model, 
we now relate it more directly to a "double solution" analogy gleaned from Couder and 
Fort [12]. These authors used the double solution ansatz to describe the behaviors of their 
"bouncer"- (or "walker"-) droplets: on an individual level, one observes particles surrounded 
by circular waves they emit through the phase-coupling with an oscillating bath, which 
provides, on a statistical level, the emergent outcome in close analogy to quantum mechanical 
behavior (like, e.g., diffraction or double-slit interference). Originally, the expression of a 
"double solution" refers to an early idea of de Broglie [ I 6] to model quantum behavior by a 
two- fold process, i.e., by the movement of a hypothetical point-like "singularity solution" of 
the Schrodinger equation, and by the evolution of the usual wavefunction that would provide 
the empirically confirmed statistical predictions. Now, it is interesting to observe that 
one can construct various forms of classical analogies to the quantum mechanical Gaussian 
dispersion [17], and one of them even can be related to the double solution idea. 

To establish correlations on a statistical level between individual uncorrelated particle 
positions x and momenta p, respectively, one considers a solution of the free Liouville equa- 
tion providing a certain phase-space distribution f (x,p,t). This distribution shows the 
emergence of correlations between x and p from an initially uncorrelated product function 
of non-spreading ( "classical" ) Gaussian position distributions as well as momentum distri- 
butions. The motivation for their introduction comes exactly from what one observes in the 
Couder experiments. In an idealized scenario, we assume that at each point x an unbiased 



emission of momentum fluctuations ttq in all possible directions takes place, thus mimicking 
(in a two-dimensional scenario) the circular waves emitted from the "particle as bouncer" . 
If we compare the typical frequency of the bouncers in the Couder experiments (i.e., roughly 
10 2 Hz) with that of an electron, for example (i.e., roughly 10 20 Hz), we see that a "contin- 
uum ansatz" is very pragmatical, particularly if we are interested in statistical averages over 
a long series of experimental runs. 

Thus, one can construct said phase-space distribution, with a being the initial x-space 
standard deviation, i.e., o"o = cr(t = 0), and 7r := mw the momentum standard deviation, 
such that 

'^-^•"{-^^H-sW- < 2i) 

Now, the above-mentioned correlations between x and p emerge when one considers the 
probability density in x-space. Integration over p provides 

p( ^ ) = / /dp = ^ exp {-# (22) 



with the standard deviation at time t given by 

a 2 = a 2 , + u 2 t 2 . (2.3) 

In other words, the distribution (2.2) with property (2.3) describes a spreading Gaussian 
which is obtained from a continuous set of classical, i.e., non-spreading, Gaussian position 
distributions of particles whose associated momentum fluctuations also have non-spreading 
Gaussian distributions. One thus obtains the exact quantum mechanical dispersion formula 
for a Gaussian, as we have obtained also previously from a different variant of our classical 
ansatz. For confirmation with respect to the latter (diffusion-based) model [4, 5], we consider 
with (2.2) the usual definition of the "osmotic" velocity field, which in this case yields 
u = —D^- = ^r. One then obtains, with bars denoting averages over fluctuations and 
positions, i.e., with x 2 = j Px 2 dx = a 
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u 2 = D 2 ( —J =D 2 I P ( -^- ) dx = ^-, and thus also u = — . (2.4) 

so that one can rewrite Eq. (2.3) in the more familiar form 

S = al (l + 5) . (2.5) 



Note also that by using the Einstein relation D = h/2m = h/4nm the norm in (2.1) thus 
becomes the invariant expression (reflecting the "exact uncertainty relation" [18]) 

1 - 1 2 
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Following from (2.5), in Grossing etal. [1, 7] we obtained for "smoothed-out" trajectories 
(i.e., averaged over a very large number of Brownian motions) a sum over a deterministic 
and a fluctuations term, respectively, for the motion in the ^-direction 



a / D 2 t 2 

x t ot{t) = vt + x(t) =vt + x{0)— =vt + x(0)\ 1+ —j-. (2.7) 

Thus one classically obtains with x (t) = x to t{t) — vt the average velocity field of a Gaussian 
wave packet as 

«**(*) = «(*) + ^r = «(*) + 1*®} $■• ( 2 - 8 ) 

For the particular choice of x (t) = o~, then, and with | — > uo for t — >■ oo, one obtains for 
large t that t>t o t(0 = ' u (^) + u o- 

Note that Eqs. (2.7) and (2.8) are derived solely from statistical physics. Still, they are 
in full accordance with quantum theory, and in particular with Bohmian trajectories [17]. 
Note also that one can rewrite Eq. (2.5) such that it appears like a linear- in-time formula 
for Brownian motion, 



x 2 = x 2 (0) + D t t } (2.9) 

where a time dependent diffusivity 



h 2 



D, = ult = ^t (2.10) 

characterizes Eq. (2.9) as ballistic diffusion. This makes it possible to simulate the dispersion 
of a Gaussian wave packet on a computer by simply employing coupled map lattices for 
classical diffusion, with the diffusivity given by Eq. (2.10). (For detailed discussions, see 
Grossing etal. (2010) [1] and Grossing etal. (2011) [ ] and the last Chapter of this paper.) 
Moreover, one can easily extend this scheme to more than one slit, like, for example, 
to explain interference effects at the double slit [7, 8]. For this, we chose similar initial 
situations as in Holland [17], i.e., electrons (represented by plane waves in the forward y— 
direction) from a source passing through "soft-edged" slits 1 and 2 in a barrier (located 
along the x-axis) and recorded at a screen. In our model, we therefore note two Gaussians 



representing the totality of the effectively "heated-up" path excitation field, one for slit 1 
and one for slit 2, whose centers have the distances +X and —X from the plane spanned 
by the source and the center of the barrier along the y-axis, respectively. 

With the total amplitude R of two coherent waves with (suitably normalized) amplitudes 
Ri = \/~Pi, and the local phases (fi, i — 1 or 2, one has, according to classical textbook 
wisdom, the averaged total intensity 



tot 



R 2 = R{ + R z 2 + 2R 1 R 2 cos<p 12 = Pi + P 2 + 2^ P 1 P 2 cos (p 12 , (2.11) 



where y9 12 is the relative phase y9 12 = y?i — tp 2 = (1*4 — k 2 ) -r. Note that y2 12 enters Eq. (2.11) 
only via the cosine function, such that, e.g., even if the total wave numbers (and thus also 
the total momenta) kj were of vastly different size, the cosine effectively makes Eq. (2.11) 
independent of said sizes, but dependent only on an angle modulo 2n. This will turn out as 
essential for our discussion further below. 

The x-components of the centroids' motions from the two alternative slits 1 and 2, 
respectively, are given by the "particle" velocity components 

v x = ±-k x , (2.12) 

m 

respectively, such that the relative group velocity of the Gaussians spreading into each other 
is given by Av x = 2v x . However, in order to calculate the phase difference ip\ 2 descriptive 
of the interference term of the intensity distribution (2.11), one must take into account the 
total momenta involved, i.e., one must also include the wave packet dispersion as described 
above. Thus, one obtains with the displacement ±x (t) = =F (X + v x t) in Eq. (2.8) the total 
relative velocity of the two Gaussians as 



Aw*, 



ot,x 



(X + V x t) 



r 2 



(2.13) 



Therefore, the total phase difference between the two possible paths 1 and 2 (i.e., through 
either slit) becomes 

I x 1 u t 

f\2 = t(mAv to t,x x) = 2mv x - - (X + v x t)x—-^-. (2.14) 

In earlier papers, Grossing (etal.) j -4], we have shown that, apart from the ordinary 
particle current J(x, t) = P(x, t)v, we are now dealing with two additional, yet opposing, 
currents 3 U = P(x, t)u, which are on average orthogonal to J, and which are the emergent 



outcome from the presence of numerous momentum fluctuations, and the corresponding 
velocities, 

We denote with u + and u_ , respectively, the two opposing tendencies of the diffusion pro- 
cess. Moreover, one can also take the averages over fluctuations and positions to obtain a 
"smoothed-out" average velocity field [2-5], i.e., 

u(M)=/pu(x, t ) d V <»<» 

In effect, the combined presence of the velocity fields u and v can be denoted as a path 
excitation field: via diffusion, the bouncer in its interaction with already existing wave-like 
excitations of the environment creates an "agitated", or "heated-up", thermal landscape, 
which can also be pictured by interacting wave configurations all along between source 
and detector of an experimental setup. Recall that our prototype of a "walking bouncer" , 
i.e., from the experiments of Couder's group, is always driven by its interactions with a 
superposition of waves emitted at the points it visited in the past. Couder etal. denote 
this superposition of in-phase waves the "path memory" of the bouncer [19]. This implies, 
however, that the bouncers at the points visited in "the present" necessarily create new wave 
configurations which will form the basis of a "path memory" in the future. In other words, 
the wave configurations of the past determine the bouncer's path in the present, whereas 
its bounces in the present co-determine the wave configurations, and thus the probability 
density distribution, at any of the possible locations it will visit in the future. Therefore, 
by calling the latter configurations a path excitation field, we point to our model's physical 
meaning of the probability density distribution: its time evolution is to be understood as the 
totality of all sub-quantum currents, which may also be described as a "heated-up" thermal 
field. 

For illustration, let us now consider a single, classical "particle" ("bouncer") following 
the propagation of a set of waves of equal amplitude Ri, each representing one of i possible 
alternatives according to our principle of path excitation, and focus on the specific role of 
the velocity fields. To describe the required details, each path i be occupied by a Gaussian 
wave packet with a "forward" momentum pj = Kki = mVj. Moreover, due to the stochastic 
process of path excitation, the latter has to be represented also by a large number iV of 
consecutive Brownian shifts as in Eq. (2.15), which on average, for stretches of free particle 

8 



propagation Vj, are orthogonal to Vj. Defining an average total velocity (with indices i = 1 
or 2 referring to the two slits, and with + and — referring to the right and the left from the 
average direction of Vj, respectively) 

v t ot,i := v, + u i+ + Uj_ (2.17) 

with the bars here (and further on, if not declared otherwise) denoting averages only over 
the spatial directions, we consider two Gaussian probability density distributions, Pi = R\ 
and P 2 = R 2 , respectively. Generally, with the total amplitude R of two coherent waves 
with (suitably normalized) amplitudes Ri = y/Pi, and the local phases cpi, i = 1 or 2, one 
has as usual that 

R = Ri cos (cot + <fx) + R 2 cos (ut + <p 2 ) , (2-18) 

which, when squared and averaged, provides the famous formula for the intensity of the 
interference pattern (2.11). Introducing an arbitrarily chosen unit vector n, one may also 
define cos (cot + cpi (x))= n • k ; (x, t), such that along with the system's evolution, the emer- 
gent outcome of the time evolution of R (x, t) as a characteristic of our path excitation field 
can analogously be assumed as 

R (x, t) = n • (r&i. (x, t) + R 2 k 2 (x, t)) . (2.19) 

Thus, with regard to the total amplitude R to t in the double-slit case, one obtains with the 
appropriate normalization that 

R 2 tot = (Xvtot,! + i? 2 ^tot, 2 ) 2 • (2.20) 

With P tot = R% ot , one can calculate the time development of the path excitation field, i.e., 
the spreading out of the total probability density in the form of the total average current 
Jtot = -PtotVtot- After a few calculational steps this provides, similarly to Grossing etal. 
(2012) [7] (but now with slightly different labellings which apply more generally), with the 
assumed orthogonality v.; • u^ = for free particle propagation (i = 1 or 2) that 



'tot 



P1V1 + P 2 v 2 + \fP~1P2 (vi + v 2 ) cos (pi2 + \fP~1P~2 (ui - u 2 ) sin ip 12 . (2.21) 



An alternative, more detailed account of Eq. (2.21) and its extension to n slits is in prepa- 

m Ri 
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ration [Fussy etal. (2013)]. Note that Eq. (2.21), upon the identification of u { = -- h x /s> 



from Eq. (2.15) and with Pi = i?|, turns out to be in perfect agreement with a comparable 
"Bohmian" derivation [17, 20]. In fact, with Vj =V5'j/m, one can rewrite (2.21) as 

3 tot = Rl^ + Rl^ (2.22) 

m m 

+ R X R 2 ( + J cos <p 12 + — (R 1 \/R 2 - R 2 Vi?i) sin <p 12 . 

\ m m J m 

The formula for the averaged particle trajectories, then, simply results from 

v tot = ^. (2.23) 

" tot 

Although we have obtained the usual quantum mechanical results, we have so far not 
used the quantum mechanical formalism in any way. However, upon employment of the 
Madelung transformation for each path j (j = 1 or 2), 



q fj = Re is / h , (2.24) 



-12 



and thus P i = R) = |^j| 2 = *^j, with the definitions (2.15) and v] := VSj/m, ipi 
(Si — S 2 )/h, and recalling the usual trigonometric identities such as cosy? = | (e llp + e~ 1,p ), 
etc., one can rewrite the total average current (2.21) immediately as 



P tot v tc 



•tot — 1 tot v tot 



(*l + * 2 )*(*l + *2)2 



i ( ife v(^i+^2) \ ^ /V v($i + ^ 



m V " (*i + ^2) / m V (^1 + ^2)* 



= -— [^*W - tf Vtf*] = —Re {tf*(-iftV)tf } , (2.25) 

2m m 

where P tot = |\&i + \I/2| 2 =: |^| 2 - The last two expressions of (2.25) are the exact well-known 
formulations of the quantum mechanical probability current, here obtained without any 
quantum mechanics, but just by a re-formulation of (2.21). In fact, it is a simple exercise 
to insert the wave functions (2.24) into (2.25) to re-obtain (2.21). 

Moreover, having thus obtained a "bridge" to the quantum mechanical formalism, one can 
now show how the Schrodinger equation not only complies with our classical ansatz, but ac- 
tually can be derived from it. Two different ways of such a corresponding derivation have al- 
ready been published by the present author (i.e., one on more general grounds [21], the other 
within a more concrete model based on nonequilibrium thermodynamics^]). Here, however, 
we just start with a result from the construction of the classical Gaussian with its dispersion 
mimicking the quantum mechanical one, i.e., more concretely, with the resulting velocity field 
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v t ot(t) from Eq. (2.8). Integrating the latter, and with £ (t) := x — vt describing the location 



of a particle in a Gaussian probability density distribution P(x,t) — ! 
we find an expression for the action function as 



27TCT 



exp 



(x-vt) 2 
2a' 2 



S 



mv tot (t)dx- 



mvx + 



■muz 



a(t) 



Edt 

1 2 

t — Et = mvx + 



mu, 



2 r 



O~0 



t- Et, 



(2.26) 



with E being the system's total energy. Note that because it generally holds that |U- = ^, 
the physics contained in (2.26) for the free particle case is essentially determined already by 
the initial conditions, i.e., where in the Gaussian the particle is initially located. Note also 
that the new term, as opposed to "ordinary" classical physics, is given by a kinetic energy, 

,1 2 



which represents a thermal fluctuation field with the kinetic temperature kT 



mu. 



m 

a(t) 



of 



what we have termed the path excitation field. Now, as with the above-mentioned Gaussian 
one has u(x,t) = —D^- = ^— = 2 °' ^0; one can rewrite the new kinetic energy term in 
Eq. (2.26) as 



mun 



k(*)J 



t 



mu 



a® 

0-Q 



t 



mu 2 



-t. 



(2.27) 



where u = u — . (This is therefore in complete accordance with the starting assumption in 



<?o 



Grossing (2008) 
that 



Moreover, upon averaging over fluctuations and space, one obtains 



P- 



muz 



a(t) 



dx 



mu}- 



flu 



(2.28) 



where the latter equation describes the identity of the averaged kinetic energy term with 
the zero-point energy as shown in Grossing (2009) [3]. 

In more general terms, i.e., independently from the particular choice of P as a Gaussian 
distribution density, with our expression for the momentum fluctuation <5p = mu as 

frVP(x,£) 



*P ( X,i)=V(*S(x, i) ): 2p(M) ■ 

we can write our additional kinetic energy term for one particle as 



^=^V(*S.V(*S) = ^(^ 



(2.29) 



(2.30) 



Thus, writing down a classical action integral for n particles, and including this new term 
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for each of them, yields (with Lagrangian L and external potential V) 



A 



L d m x dt 
P 



dS 
~dt 



v^ i „ n „ n A i fhViP\ 



V 



d m xdt, 



(2.31) 



where P = P(xi,x 2 , . . . ,x n , £). Using again the Madelung transformation (2.24) where 
R = \/P, one has 

(2.32) 



V^ 



V> 



d m a;dt 



V^ 



^ 



2 P 



/j 



and one can rewrite (2.31) as 

A= [ Ldt = I <Y".rdt 



2 fdS 



'^U 



\ n 

+ ")+£ 

7 i=i 



2m, 



|V^| : 



(2.33) 



Thus, with the identity (V'Pff = —^(ip*ip — ip*ip) : one obtains the familiar Lagrange density 

n %2 



\h n 2 

1=1 



(2.34) 



from which by the usual procedures one arrives at the n-particle Schrodinger equation 



ift 



dip 
~dt 



h 



E^+^- 



i=l 



(2.35) 



Note also that from (2.31) one obtains upon variation in P the modified Hamilton- Jacobi 
equation familiar from the de Broglie-Bohm interpretation, i.e., 



as 

~dt 



+e 

i=\ 



2rrii 



+ y(xi,x 2 ,...,x n ,t) + [/(xi,x 2 ,...,x n ,£) = 0, 



(2.36) 



where U is known as the "quantum potential" 



£7(xi,x 2 ,...,x n 



>') = £ 



8=1 



fi 2 

4m,- 



1 /VjP 

2 V P 



VfP 
P 



£ 



»=i 



fi 2 V 2 P 

2m,- P 



(2.37) 



Moreover, with the definitions 



Spi h v,p , . iViP v,p 

u, := = and k,,, := = 

m,- 2m, P 2 P P 



one can rewrite U as 

n 

^ = E 



»=i 



mjUj ■ Uj h 



(Vi • Ui) 



E 



»=i 



2m, 



A-"-ut ' *-m V, ■ k U j 



(2.38) 



(2.39) 
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However, as was already detailed in Grossing (2009) [3], for the energy balance kT = hu> 
referring to the vacuum's acting as the particle's "thermostat" , Uj can also be written as 

Ui = ^—VtQ , (2.40) 

which thus explicitly shows its dependence on the spatial behavior of the heat flow Q = 
kT In P. Insertion of (2.40) into (2.39) then provides a thermodynamic formulation of the 
quantum potential as 

~1 fViQ\ 2 V^Q 



h 2 



Am; 



hwi ) huji 



(2.41) 



3. "SYSTEMIC" NONLOCALITY IN DOUBLE SLIT INTERFERENCE 

In a recent paper, Tollaksen etal. [ '] renew the discussion on interference at a double 
slit "from a single particle perspective", asking the following question: If a particle goes 
through one slit, how does it "know" whether the second slit is open or closed? We shall 
here first recapitulate the arguments providing these authors' answer and later provide our 
own arguments and answer. Of course, the question is about the phase information and how 
it affects the particle. We know from quantum mechanics that phases cannot be observed on 
a local basis and that a common overall phase has no observational meaning. Assuming that 
two Gaussian wave functions, \l/i and \&2, describe the probability amplitudes for particles 
emerging from slits 1 or 2, respectively, which are separated by a distance D, the total wave 
function for the particle exiting the double slit may be written as 

# = e iai # 1 + e ia2 ^ 2 , (3.1) 

but since a common overall phase is insignificant, one multiplies (3.1) with e~ 1Ql and writes 
the total wave function as 

^ = ^i + e^^ 2 , (3.2) 

where ip := «2 — «i is the physically significant relative phase. Tollaksen etal. now ask 
where the relative phase appears in the form of (deterministic) observables that describe 
interference. For the following recapitulation of their argument, it must be stressed that the 
authors assume the two wave functions, \l/i and ^ 2 , to be initially non-overlapping. When 
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looking at the expectation value (in the one-dimensional case, for simplicity) 

x = f %x^ 9 dx = J (*; + e -1 ^) x (*! + e 1 ^) dx (3.3) 

= /(*; rfl + ^ 2 ,d. + /*;^ 2dl + c,., 

one sees that it is independent of p because of the vanishing of the last term due to the 
assumed non-overlapping of \?i and ^2- Similarly, this also holds for all moments of x, 
and for all moments of p as well. In particular, one has for the expectation value of the 
momentum 

%ih—* v dx (3.4) 

= Re J (^* 1 ih—^ 1 + ^* 2 ih—^ 2 j dx + Re f ^>]ih—e^ 2 dx + c.c, 

where the ^-dependent term vanishes identically, because -§^2 = for ^ 2 = 0. So, 
again, where does the relative phase appear? The answer of Tollaksen etal. is given by 
a "shift operator" that shifts the location of, say \&i, over the distance D to its new location 
coinciding with that of \&2- The expectation value of the shift operator is thus given by 



e - i£ ? = / Ufle-'f^x + %e' iE ^^ 2 j dx (3.5) 

+ f^* 1 e- iE ^^ 2 e i ' p dx+ /"^e^e-^^idx, 
where all but the last term vanish identically, thus providing (with the correct normalization) 



e 1 h = e ltp /2 and e 1 » + e 1 r = cos (p. (3.6) 

In order to physically interpret the shift operator, the authors now shift to the Heisenberg 

2 

picture, thereby providing with a Hamiltonian H = ^— + V (x) the time evolution of the 
operator as 



d -pd i 
— e h = — 
dt h 



,2 



2m 



\D _ iE R (V{x)-V{x + D) 



e-V w v _ (3 _ 7) 



h \ D 

With its dependence on the distance D between the two slits, Eq. (3.7) is a description 
of dynamical nonlocality: it is thus shown how a particle can "know" about the presence 
of the other slit. Tollaksen etal. maintain that it is possible to understand this dynamical 
nonlocality only by employing the Heisenberg picture. However, we shall now show that such 
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an understanding is possible also within the Schrodinger picture, and even more intuitively 
accessible, too. For, there is one assumption in the foregoing analysis that is not guaranteed 
to hold in general, i.e., the non-overlapping of the wave functions \I/i and ^ 2 . On the contrary, 
we now shall assume that the two Gaussians representing the probability amplitudes for the 
particle immediately after passing one of the two slits do not have any artificial cut-off, but 
actually extend across the whole slit system, even with only very small (and practically 
often negligible) amplitudes in the regions further away from the slit proper. (We shall give 
arguments for this assumption further below.) In other words, we now ask: what if \Pi and 
^2 do overlap, even if only by a very small amount? To answer this question, we consider 
the expectation value of the momentum and obtain from Eq. (3.4) that the terms involving 
the relative phase if provide with \l/j = Rje lifj and tpj = -j- 

p = R X R 2 (VSi + VS 2 ) cosip + h (R 1 VR 2 - R2VR1) sin <p. (3.8) 

First of all one notes upon comparison of (3.8) with Eqs. (2.21)-(2.23) the exact correspon- 
dence of p with our classically obtained expression for the interference terms of the emerging 
current, or of the expression for mv M) respectively. Moreover, although the product R±R 2 
is in fact negligibly small for regions where only a long tail of one Gaussian overlaps with 
another Gaussian (i.e., such that the non-overlapping assumption would be largely justified), 
nevertheless the second term in (3.8) can be very large despite the smallness of R\ or R 2 . 
It is this latter part which is responsible for the genuinely "quantum mechanical" nature of 
the average momentum, i.e., for its nonlocal nature. This is formally similar in the Bohmian 
picture, but here given a more direct physical meaning in that this last term refers to a dif- 
ference in diffusive currents as explicitly formulated in the last term of Eq. (2.21). Because 
of the mixing of diffusion currents from both channels, we call this dominant term in J tot 
the "entangling current" [ 3], J e . In fact, inserting Eq. (2.40) for a one-particle system into 
the latter expression, one obtains 

t /73-73-V (Qi — Q2) . (Q n N 

J e = \JP1P2 o sin<pi2. (3.9) 

zoom 

In other words, the entangling current is essentially a measure for the gradient of the vac- 
uum's "heat" fluxes Qi associated with the channels i. 

Finally, there are substantial arguments against the non-overlapping scenario in Tollak- 
sen et al. Firstly, experiments by Rauch et al. have shown that in interferometry interference 
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does not only happen when the main bulks of the Gaussians overlap, but also when a Gaus- 
sian interferes with the off-bulk plane-wave components of the other wave function as well 
1 24]. On a theoretical side, we have repeatedly stressed that the diffusion processes em- 
ployed in our model must be described by nonlocal diffusion wave fields [25, 26] which thus 
require small but non-zero amplitudes across the whole experimental setup. Moreover, and 
more specifically, we have shown [ , ] that quantum propagation can be identified with 
sub-quantum anomalous (i.e., "ballistic") diffusion which is characterized by infinite mean 
displacements x = oo despite the finite drifts x 2 = u 2 t 2 and u < c. In sum, these arguments 
speak in favor of using small, but non-zero amplitudes from the Gaussian of the "other slit" 
which can interfere with the Gaussian at the particle's location. 

Let us now see how we can explain dynamical nonlocality within our sub-quantum ap- 
proach. For example, we can ask the following question in our context: what if we start 
with one slit only, and when the particle should pass it we open the second slit? Let us 
assume for the time being, and without restriction of generality, that the ^-component of 
the velocity v x is zero. Then, according to (2.14), upon opening the second slit (with the 
same a), we obtain a term proportional to the distance 2X between the two slits. Thus, 
there will be a shift in momentum on the particle passing the first slit given by 

^4v^±±V( Sl - S2 H%i, (3,0) 

where one effectively uses the "modular momentum" p mo d = p mod^ = p — lnit^ because 
an added or subtracted phase difference y>\i = 2nir does not change anything. In other 
words, by splitting X into a component X n providing ip±2 = Iwn on the one hand, and the 
modular remainder AX on the other hand providing ip\i = p mo dx/h, one rewrites (2.14) 
with X := X n + AX (and with v x = for simplicity) as 

l u 2 t h t 

¥?i2 = ~{X n + AX)x- — =: -2nn - AXx- ^ . (3.11) 

D ctq 2m cHctq 

Therefore, one can further on substitute Ap by Aj9 mot j, and one obtains a momentum 

shift (with the sign depending on whether the right or left of the two slits is opened as the 

second one) 

h AX h 2 t A D 2 t A & . 

A Pmod = ±-V Vl2 = ±^^ 5-2 = ±mAX^^ = ±mAX- . 3.12 

2 2 2m cHctq a G t o 

It follows that, due to the "vacuum pressure" stemming from the opened second slit, there is 

an emergent nonlocal "force" which does not derive from a potential but from the impinging 
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of the second slit's sub-quantum diffusive momenta on the particle at the first slit. As there 
is no additional force on the particle, we do not speak of a "dynamical" nonlocality, but of 
a "systemic" one instead. For more details, see Grossing etal. (2013) [ ]. 

4. QUANTUM MECHANICAL RESULTS OBTAINED FROM CLASSICAL COM- 
PUTER SIMULATIONS 

For our classical simulations of quantum phenomena we make use of the physics of ballistic 
diffusion as given in Eq. (2.9), with the time-dependent diffusivity D t from Eq. (2.10). We 
use an explicit finite difference forward scheme [28], 

OP 1 

^^0PM + i]-PM), (4.i) 

with 1-dimensional cells. For our case of D t being independent of x, the complete equation 
after reordering reads as 



D[t + 1]M 

~Ax~ 2 



P[x, t + 1] - P[x, t] + L 7 2 — i P l x + M] ~ 2P[x, t] + P[x - 1, t}} (4.3) 



with space x and time t, and initial Gaussian distribution P(x, 0) with standard deviation 
co- 
in the following, some exemplary figures with one spatial and one time dimension are 
shown to demonstrate some of the present applications of our model as well as of the novel 
simulation protocol that comes along with it. The first three figures exhibit double-slit 
scenarios, with time evolving from bottom to top, and the remaining three figures show 
n-slit scenarios with time running from left to right. Except for the last figure, the coloring 
for intensity distributions is, for increasing intensity, from white through yellow and orange, 
and average trajectories are displayed in red. 

Turning now to Fig. 4.1, we begin with the upper part, i.e., a quantum mechanical interfer- 
ence pattern obtained from a computer simulation employing classical physics only. Shown 
are the probability density distributions for two wave packets emerging from the Gaussian 
slits with opposite velocities, superimposed by flow lines which coincide with Bohmian tra- 
jectories. Note, however, that the latter appear only as the result of the averaging, whereas 
the sub-quantum model on which the figure is based involves wiggly, stochastic motions. 
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For the smoothed-out trajectories, however, a "no crossing" rule applies with respect to the 
central symmetry axis, just like in the Bohmian picture. The figure's lower part shows the 
"entangling current" of Eq. (3.9) for the scenario above, which is in our model the decisive 
part of the probability density current that is responsible for the genuinely "quantum" ef- 
fects. Due to the sinusoid nature, colors are used that display positive (red) and negative 
(blue) values. Note that after the time of maximal overlap between the two wave-packets, 
the order of maxima and minima is reversed, which results from the respective differences 
in the exchanges of "heat" according to Eq. (3.9). 

Fig. 4. 2 shows a similar scenario as Fig.4.1, albeit with zero velocity components in the 
transverse (i.e., x— ) direction. However, whereas the situation of Fig.4.1 is characterized by 
small dispersion of the Gaussians, here the dispersion is significantly higher. The effects of 
the "no crossing" rule are clearly visible. Note also the "kinks" of trajectories moving from 
the center-oriented side of one relative maximum to cross over to join more central (relative) 
maxima. In our model, a detailed micro-causal account of the corresponding kinematics can 
be given. As can be seen from the entangling current in the lower part of the figure, its 
extrema coincide with the areas where the kinks appear, pointing at a rapid crossing-over 
due to the effects of the diffusive processes involved. 

Fig. 4. 3 shows almost the same scenario as Fig. 4. 2, the only difference being that in the 
left channel a phase shift of Aip = | has been added, essentially leading to a representation 
of the scalar Aharonov-Bohm effect. Note that the minima and maxima of the intensity 
distribution are shifted accordingly. This comes along with a different behavior of the 
entangling current. Whereas in Fig.4.2 the entangling current is asymmetrical with respect 
to the central symmetry line, it now is symmetrical. Effectively, this means that due to the 
diffusion processes providing some "vacuum pressure" , the intensity maxima are "pushed to 
the side", providing an asymmetrical distribution of the interference fringes. 

Fig. 4. 4 displays two examples which extend the double slit simulations to more than 
two slits. Whereas the upper picture shows a 3-slit system with intensity distributions and 
averaged trajectories, the lower one is a detail of a "Talbot carpet", i.e., a system with a 
large number n of slits, albeit exhibiting only four of them. It clearly demonstrates a result 
that is well-known from comparable Bohmian calculations, i.e., that for the case of large n, 
non-spreading "cells" emerge which contain the trajectories and keep them there [. ]. 

To continue with ra-slit systems, Fig. 4. 5 shows, in the upper part, a type of Talbot carpet 
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(without trajectories) with gradually diminishing initial intensity distributions, from top to 
bottom. The lower figure shows the same except that the initial intensity distributions have 
been chosen to have random values. This is an example of a system whose complexity is 
considerable, despite the fact that the simulation is still very fast and simple, thus illustrating 
also an advantage over purely analytical tools. Eventually, this type of simulations might 
even apply at levels of complexity which would exceed any traditional analytical methods 
of quantum mechanics. 

Finally, Fig. 4. 6 shows two examples of intensity and trajectory distributions of partic- 
ularly "weighted" 7-slit systems. In the upper figure, the high relative intensities at the 
extreme locations force the trajectories of the remaining slits to become "squeezed" into a 
narrow canal. The lower figure instead shows the opposite effect of a centrally positioned 
relatively high intensity "sweeping" away the other slits' trajectories. The two examples 
again illustrate the effect of the vacuum's "pressure" in our sub-quantum model. 
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Figure 4.1. Top: Quantum mechanical interference pattern, with averaged trajectories, obtained 
from a computer simulation employing classical physics only. Bottom: "Entangling current" of Eq. 
(3.9) for the scenario above. 
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Figure 4.2. Similar scenario as in Fig. 4.1, albeit with zero velocity components in the transverse 
(i.e., x—) direction, and with large dispersion instead of a small one. As can be seen from the 
entangling current in the lower part of the figure, its extrema coincide with the areas where the 
trajectories' "kinks" appear, pointing at a rapid crossing-over due to the effects of the diffusive 
processes involved. 
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Figure 4.3. Same as Fig. 4. 2, except that in the left channel a phase shift of A92 = | has been 
added, thus demonstrating an example of the scalar Aharonov-Bohm effect. Due to the entangling 
current, or to the diffusion processes, respectively, the intensity maxima of Fig. 4. 2 are "pushed to 
the side" , providing the asymmetrical distribution of the interference fringes. 
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Figure 4.4. Top: A 3-slit system with intensity distributions and averaged trajectories. Bottom: 
Detail of a system with a large number n of slits (i.e., a so-called Talbot carpet), explicitly showing 
only four of them. Note that non-spreading "cells" emerge which contain the trajectories and keep 
them there. 



25 



Figure 4.5. Top: Classically simulated interference pattern for a system of 9 slits with gradu- 
ally diminishing initial intensity distributions. Bottom: Same, except that the initial intensity 
distributions have been chosen to have random values. 
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Figure 4.6. Two examples of intensity and trajectory distributions of particularly "weighted" 7-slit 
systems. Top: Quantum squeezer. Bottom: Quantum sweeper. 
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